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Abstract 

Here we consider a class of 2 (g) 2 d chessboard density matrices starting with three- 
qubit ones which have positive partial transposes with respect to all subsystems. To 
investigate the entanglement of these density matrices, we use the entanglement witness 
approach. For constructing entanglement witnesses (EWs) detecting these density ma- 
trices, we attempt to convert the problem to an exact convex optimization problem. To 
this aim, we map the convex set of separable states into a convex region, named feasible 
region, and consider cases that the exact geometrical shape of feasible region can be 
obtained. In this way, various linear and non-linear EWs are constructed. The optimal- 
ity and decomposability of some of introduced EWs are also considered. Furthermore, 
the detection of the density matrices by introduced EWs are discussed analytically and 
numerically. 

Keywords: chessboard density matrices, optimal non-linear entanglement 
witnesses, convex optimization 
PACs Index: 03.65.Ud 



Three-qubit EWs 

1 Introduction 



3 



Bound entangled states, states with positive partial transposes with respect to all subsystems, 
are of great importance in quantum information processes ^ El El IH E] • One class of bound 
entangled states is the three-qubit chessboard states considered in [6] where the authors have 
used a separability criterion due to P. Horodecki to show the boundness of such states. The 
boundness of these states for some range of parameters are also investigated in [7J using 
entanglement witnesses (EWs) and in [8] from the perspective of convex optimization. Another 
class of chessboard states has been discussed in |9] again by using entanglement witnesses 
(EWs). The EWs are of special interest since it has been proved that for any entangled state 
there exists at least one EW detecting it. The EWs are Hermitian operators which have 
non-negative expectation values over all separable states while they have negative expectation 
values over, that is they are able to detect, some entangled states [TOl ITT] . 

In this paper, we consider a generalized form of the above chessboard states initially for 
2 C?) 2 ® 2 case, then extend them for 2 ® 2 ® d case and use EWs approach to analyze their 
entanglement. For constructing the relevant EWs, we attempt to convert the problem to 
an exact convex optimization problem. This method are general and one can apply it for 
multi-qubits in a similar way. All of witnesses constructing in this way are valid with some 
changes in notation. As the dimension of problem increases the number and categories of 
EW's increases but the procedures are same in general. Convex optimization techniques have 
been widely used in quantum information problems recently [121 [El El [ISl [161 113 HI Ell 1201 
EH |22l El]. In references [Ml ESI [26l [271 [2H] the problem of constructing EWs was converted 
to a linear programming problem, a special case of convex optimization problem, exactly or 
approximately. To this aim, the convex set of separable states was mapped into a convex 
region, named feasible region (FR). The FR may be a polygon by itself or it may not. When 
FR was not a polygon, it was approximated by a polygon. In this way, the problem was 
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converted to a linear programming problem whose linear constraints came from the exact or 
approximated boundary surfaces of FR. 

Here we consider the cases that the geometrical shape of FR can be obtained exactly and 
hence convert the problem to an exact convex optimization problem. Any hyper-plane tangent 
to the FR corresponds to a linear EW. According to the geometrical shape of FR, we can 
construct non-linear EWs or can not. It is shown that when the geometrical shape of FR is 
a polygon, all EWs are linear; otherwise it is possible to construct non-linear EWs. In the 
previous works where a non-polygonal FR was approximated by a polygonal one, the number 
of obtained linear EWs was not sufficient for constructing non-linear EWs. However, in the 
present work where we consider the exact geometrical shape of a non-polygonal FR, any hyper- 
plane tangent to the surface of FR is a linear EW. Therefore, there exist innumerable linear 
EWs which is enough for constructing a non-linear EW as the envelop of linear functionals 
arising from them. By construction, a non-linear EW plays the role of innumerable linear EWs 
as a whole and hence it may detect bound entangled states. Our approach is typical and can 
be applied in all cases where the exact geometrical shape of FR is known. 

The paper is organized as follows. In Section 2, we review the basic notions and definitions 
of EWs relevant to our study and describe our approach of constructing EWs. Then we present 
a generalized form of a class of three-qubit density matrices of [6]. In Section 3, we consider 
the construction of linear and non-linear EWs that can detect the mentioned density matrices. 
Section 4 is devoted to an analysis of optimality of introduced EWs. It is proved that some 
of the EWs are optimal. In Section 5, we consider the detection of mentioned density matrix 
by introduced EWs analytically and numerically. Section 6 is devoted to the comparison of 
our results with other works. In section 7 we extend all these methods to 2 (g) 2 ® (i case and 
we see that all these methods are general and one can apply them for multipartite chessboard 
density matrices. This extension neither change the structure of PPT's conditions nor the 
EW's structures. In section 8 numerical analysis for detection ability of introduced EW's for 
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2® 2® 2 and 2 ® 2 3 chessboard density matrices are discussed. 
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2 Preliminaries 

2.1 A class of three-qubit density matrices with positive partial transposes 

Here we consider a generalized form of a class of three-qubit density matrices presented in [6] 
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(2.1) 



where a, 6, c, d are non-negative parameters, < < 1 for i = 1, 2, 3, 4 and n = { a + h + c + 
d + \ + \ + \ + \ )- It is easy to see that this density matrix has positive partial transposes 
with respect to all subsystems, i.e., it is a PPT state. The density matrix of [6] is a special 
case of p where = 0, = 1, = = = 0, and a = 1. We want to show that for some 
values of the parameters, p is a PPT entangled state. To this aim, we will construct various 
linear and non-linear non-decomposable EWs that are able to detect it. 
Written in the Pauli matrices basis, p has the form 

p = III 1 1 + r^^^ajl + r^^jaj + r^^JIa, + r^^^^a.aj + r^^^aja, 

+rg^^Iazaz + Vyy^o^o^o^ + r ^^^o ^^a + r^^^o^o^Oy + r^^^o^OyO^ (2.2) 

~t~^2ii ^y^a^^a^ ~^ ^2i2^y^x^y ~t~ ^22i^y^y^x ~^ '^222^y^y^y\ 

where the coefficients r.^^ are given in appendix B. We will try to construct our non-decomposable 
EWs by using Pauli group operators appearing in the p. But before this, let us review the 
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2.2 Entanglement witnesses 

Let us first recall the definition of entanglement and separability [29]. By definition, an n- 
partite quantum mixed state p G B{T-C) (the Hilbert space of bounded operators acting on the 
Hilbert space H = Hd^ ... Ti-d,,) is called fully separable if it can be written as a convex 
combination of pure product states, that is 

P = Y:pM'^){4^\ ® l«f ® - ® (2.3) 

i 

where are arbitrary but normalized vectors lying in the 'Hdj, and Pi > with X^jP* ~ 

1. Otherwise, p is called entangled. Throughout the paper, by separability we mean fully 
separability. 

An entanglement witness (EW) W is a Hermitian operator which has non-negative expec- 
tation value over all separable states ps and its expectation value over, at least, one entangled 
state pe is negative. The existence of an EW for any entangled state is a direct consequence 
of Hahn-Banach theorem [30] and the fact that the subspace of separable density operators is 
convex and closed. 

Based on the notion of partial transpose map, the EWs are classified into two classes: 
decomposable (d-EW) and non-decomposable (nd-EW). An EW W is called decomposable if 
there exist positive operators V, Qk such that 

W = r+^ Q]^ (2.4) 

where M := {1, 2, 3, n} and Tk denotes the partial transpose with respect to partite K G Af 
and it is non-decomposable if it can not be written in this form [3T|. Clearly, d-EWs can not 
detect bound entangled states (entangled states with positive partial transpose (PPT) with 
respect to all subsystems) whereas there are some bound entangled states which can be detected 
by an nd-EW. 
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A non-linear EW associated to an entangled density matrix g is simply a non-linear func- 
tional of g such that it is non-negative valued over all separable states, but has negative value 
over the density matrix g. A non-linear EW can be viewed as the envelop of a set of linear 
functionals Tr{Wg) that arise from corresponding linear EWs W. 

Usually one is interested in finding EWs W which detect entangled states in an optimal 
way. An EW W is called an optimal EW if there exists no other EW which detects more 
entangled states than W. It is shown that the necessary and sufficient condition for optimality 
of an EW W is that there exist no positive operator V and e > such that W = W — eV he 
again an EW [32]. 

2.3 Constructing of EWs via exact convex optimization 

Let us consider a set of given Hermitian operators Q-. By using these operators, we will 
attempt to construct various linear and non-linear EWs. To this aim, for any separable state 
Ps we introduce the maps 



which map the convex set of separable states into a convex region named the feasible region 
(PR). Any hyper-plane tangent to the FR corresponds to a linear EW, since such hyper-planes 
separate the FR from entangled states. Hence, we need to determine the geometrical shape 
of FR. In general, determining the geometrical shape of FR is a difficult task. However, one 
may choose the Hermitian operators Q- such that the exact geometrical shape of FR can be 
obtained rather simply. By such a choice, when the FR is a polygon, its surface corresponds 
to linear EWs which are linear combinations of the operators Q- ; otherwise, linear EWs come 
from any hyper-plane tangent to the surface of FR. When the FR is not a polygon, besides 
the linear EWs it is possible to obtain non-linear EWs for the given density matrix. 

To obtaine the geometrical shape of FR, we note that every separable mixed state ps can 
be written as a convex combination of pure product states, so the subspace of separable states 
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S can be considered as a convex hull of the set of all pure product states V. Thus first we 
specify the geometrical shape of a region obtained from mapping of V under the Pj's. If the 
resulted region is convex by itself, we get the FR, otherwise we have to take the convex hull 
of that region as FR. 

In this paper, the operators are chosen as linear combinations of Hermitian operators 
in the Pauli group Qn, a group consisting of tensor products of the identity I2 and the usual 
Pauh matrices a^, Cy and together with an overall phase ±1 or ±i [231 EH 155] . 



3 A class of three-qubit EWs 

In this section,we want to introduce various nd-EWs for the density matrix p of (12.11) . To 
simplify the analysis, let us classify these EWs according to the shape of relevant FRs: polyg- 
onal, conical, cylindrical and spherical. Hereafter, we will use the following notation for the 
three-qubit Pauli group operators 

Oijk = cTj ® (Tj (g) (Tfc, i, j, k = 0,1, 2, 3, (3.6) 

where ctq, cti, (72 and stand for the 2x2 identity matrix I2 and single qubit Pauli matrices 
cr^, ay, (T2 respectively. Let us begin with polygonal case. 

3.1 EWs with polygonal FR 

Let us consider the following operators 

qT = O333, qT = Oni + {-iyoi22, qT = O212 + (-i)*+'022i, ^ = 0, 1 

and try to construct nd-EWs from them for detecting p. To this end, we define the maps 

P, =Tr(Q;»|/3)|7)(«|(/5|(7l), J = 1,2,3 
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for any pure product state |q;)|/3)|7). In this case, the FR is a polygon which its boundary 

planes are as follows: 

{-iy.p, + {-iy.p, + {-i)hp^^i , (j\,j„j3)e{o,ir (3.7) 

(for a proof, see appendix A). These planes can be rewritten as 



min Tr([III - (-ly^Ql" - {-ly^Ql - {-iy^Ql"]\am\^){am{^\) = 

I«>l/3>|7> 



It is seen that the operators in the bracket have non-negative expectation values over all pure 
product states, hence they give rise to the following hnear EWs 

'Klv, = ^m(-l)'^0333+(-l)'^Olll + (-l)'^Ol22 + (-l)'^0212 + (-l)'^ +'3+^4+10221, (3.8) 

where e {0, 1}^. Besides the above 16 EWs, we can construct other 16 EWs by 

using the fact that local unitary operators take an EW to another EW. For this purpose, we 
act the phase-shift gate 



^ 0^ 



locally on the first qubit which takes — > a^, — > ^-iid — under conjuga- 
tion, and get 

Xlv. = ^^mZ^.,,Wn = III + (-1)^^0333 

+ (-l)*^Om + (-l)^30i22 + (-1)*4+'0212 + (-1)^2+^3+^40221. 

We could replace with the operator azCr^I or any cyclic permutation of it, but since these 
lead to d-EWs we do not consider such cases here. 

In this way, we have constructed 32 hnear EWs with polygonal FR. 

3.2 EWs with conical FR 

For this case, we consider the following Hermitian operators 

QT-Ok'j'i', qT ^Om + i-iyokji, qT ^Oikj + {-iyOjik, i = o,i. 
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where k'j'l' is one of the triples 333, 330, 303, 033, and kjl is one of the triples 122, 212, 221. 
Now we try to determine the exact shape of the FR. The FR is a cone given by 

(1 ±Pi)^ = P| + P| (3.10) 

(for a proof, see appendix A), where 

P=Tr(Q;»|/5)|7)(«|(/3|(7|), J = 1,2,3. 

We assert that any plane tangent to the FR corresponds to an EW. To show this, we maximize 
the function 

f{P,,P2,P3)=A,P,+A,P2 + A,P3 (3.11) 

where A- are real parameters, under the constraint f l3.10p . This is a convex optimization 
problem since the function and its constraint are both convex functions. Using the Lagrange 
multiplier method shows that this maximum is ±^4^ provided that = + A^. It is easy 
to see that the plane A^Pi + A^P2 + A.^Ps = ±A-^ is tangent to the surface (13.10p at the point 
(— ± 1,A,^,A.^). This plane can be rewritten as 

min Tr{[AJII±A,Q';" ± (A,Q^° + A,Q'';)]\am\^){am{^\) = 0. 

Thus the operator 

W^" = AJII ± A,gf° ± (A,Qj^° + A3Q^°) 

has non-negative expectation value over all pure product states, hence it can be a linear EW. 
By defining cosip = ^ and sinip = W^" is rewritten as 

''''''W^S,^± = III ± Ou'fv + cos^ (Oni + {-irOu.i) + sinV^ (0,^, + (-1)^0,,^) . (3.12) 

where z = 0, 1. Now we obtain non-linear functionals of p, hence non-linear EWs, by optimizing 
Tr[{^'^'^'w'^°nj^p] with appropriate choice of the parameter ■?/' as a functional of p. We note 
that 

Trt''''wt;,^)p\ = 1 ± r^,^,„ + cosz^(r,,, + (-1)V,^J + sin^(r,,^. + (-1)V^,J. 
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By defining 

ill ' V ' I,. ,j 

COS rj — 



^(r,,, + (-l)V,J2 + (^^^^ + (_l)V^j. 
Tr[{^'^'^'w'^°l^^J^p\ can be rewritten as 



The trace take its minimum for tjj — t] = tt: 



'''''F^;,i±iP) - ^^Tr[e'^''Ki.±)p] = 1 ± r,,,, - V(r,,, + (-l)X,J2+(v. + (-l)V,J2. 

(3.13) 

These arc the required non-hnear functionals, hence non-hnear EWs, associated with p. It is 
seen that the number of such non-hnear EWs is 48. 

We can obtain other 48 hnear EWs from " ^ ' ^kji i± conjugating them with MIL This 
gives further 48 non-hnear EWs of conical case as foUows 



(3.14) 



here kjl is one of the triples 211, 121 and 112. 

In this way, we have constructed 96 non-linear EWs with conical PR. 

3.3 EWs with cylindrical FR 

The second type of non-linear EWs for p can be derived by considering the following operators 

QT-Ok'ri', QT ^O^n + {-lyOkjl, QT = Oikj + {-iy^'Ojik, 1 = 0,1, 

where k'fl' is one of the triples 300, 030, 003, and kjl is one of the triples 122, 212, 221. It 
can be shown that the FR has the cylindrical shape 

+ (P2 + P3)' = 1, (3.15) 
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the maximum of the function (13.111) under the constraint f l3.15p is A?^^ A?^ provided that 
= A^ and this leads to the hnear EWs 

J ' 1 2 

(3.16) 



where cosip = / i^A^ + and i^ji^ = 0,1. Similar arguments as above shows that 
^'^'^'W^^L.^ j gives rise to non-linear EWs for p as follows 

■''12 

The number of these non-linear EWs is 36. We obtain other 36 non-linear EWs of this type 
by conjugating ^'^'^'Wkji^i^i^ with Mil as follows 

= minTr[M//('='^''V5,,^,jMtjJp] 
= 1 - ^Jrl^,, + (r,,, + (-l)^^r,^, + (-l)^^r,,, + (-l)H+^.+ir ^ 
here kjl is one of the triples 211, 121 and 112 and i^ii^ = 0, 1. 

In this way, we have constructed 72 non-linear EWs with cylindrical PR. 

3.4 EWs with spherical FR 

The third type of non-linear EWs for p follows from the operators 

where k'j'V is one of the triples 300, 030, 003, and kjl is one of the triples 122, 212, 221. In 
this case, the FR is of spherical shape 

Pl + P^ + Pl = l, (3.19) 

the maximum of the function (13.111) under the constraint (13.19^ is ^J A^^ + A^ -\- and this 
leads to the linear EW 

^'^'^'^Tj^ = III + {smr]cosC)Ok'fi' + sinr/ sin ((Om + {-lyOkji) + cosr]{Oikj + {-lyOjik), 

(3.20) 
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where 

sin 77 cos C= — ^ , sin 77 sin C= — ^ , cos 77 = — 

^A'i + Ai + Af ^Ai + Ai + A{ ^A{ + + Ai 

The 18 non-hnear EWs which correspond to ^'^'^'W^^^ ^ is 

''^■''>S;.(P) = min Tr\{^'''''W%,M = 1 - ^^f,^.,, + (^n + + (r,,, + (-1)V,,J2. 

(3.21) 

We obtain other 18 non-hnear EWs of this type by conjugating "^W^^ with Mil as follows 

(3.22) 



= 1 - + (^... + (-1)^^.^' + {r,,, + (-1)X,J^. 

here /cj7 is one of the triples 211, 121 and 112. 

In this way, we have constructed 36 non-hnear EWs with spherical PR. 



4 Optimality of the EWs 

In this section we discuss the optimality of EWs introduced so far. Let us recall that if there 
exist e > and a positive operator V such that W = W — eV he again an EW, the EW 
W is not optimal, otherwise it is. Every positive operator can be expressed as a sum of pure 
projection operators with positive coefficients, i.e., V = AJ'i/'.)(z/'J with all > 0, so we 
can take V as pure projection operator V — If W is to be an EW, then l'^) must be 

orthogonal to all pure product states that the expectation value of W over them is zero. The 
eigenstates of each three-qubit Pauli group operator can be chosen as pure product states, half 
with eigenvalue -|-1 and the other half with eigenvalue -1. In EWs introduced so far, there 
exists no pair of locally commuting Pauli group operators, so the expectation value of such 
pauli group operators vanishes over the pure product eigenstates of one of them. 

Regarding the above facts, now we are ready to discuss the optimahty of introduced EWS. 



Three-qubit EWs 14 
4.1 Optimality of EWs with polygonal FR 

Let us begin with EWs of (13. 8p . We discuss two cases = and = 1 separately. For the 
case i-^ = 0, note that as eigenstates of the operator 0-^0-^0-^ with eigenvalue +1 we can take 
the pure product states 

\z;+)\z;+)\z;+), \z- +)\z- -)\z; -), \z; -)\z; +)\z; -), \z; -)\z- -)\z; +) , (4.23) 

and as eigenstates with eigenvalue -1 we can take the following ones 

\z-+)\z-+)\z--), \z;+)\z;-)\z;+), \z; -)\z; +)\z- +) , \z; -)\z; -)\z- -) . (4.24) 



The EWs W^oijjgi^ have zero expectation values over the states of (14.241) . so if there exists a 
pure projection operator that can be subtracted from EWs VFQ^j°j^j^, the state ought 

to be of the form 

1^) = a+++\z; +)\z; +) + a^__\z; +)\z; -)\z; -) 

(4.25) 

+a_+_\z; -)\z; +)\z; -) + a__Jz; -)\z; -)\z; +). 

Expectation values of W^°^j^ over pure product eigenstates of the operator (TxCTxCrx with eigen- 
value -1 are zero, so should be orthogonal to these eigenstates. Applying the orthogonality 
constraints gives the following equations 

(x; +\{x; +\{x; -Wtp) = - - + a__+) = 0, 

{x; +\{x; -\{x; = - + - a__+) = 0, 

{x; -\{x; +\{x; = + a^__ - a_^_ - a__+) = 0, 

{x; -\{x; -\{x; = + a+__ + a_^_ + a__+) = 0. 

The solution of this system of four linear equations is a^_^_^ = = a_^_ = a ^ = 0. Thus 

Itp) = 0, that is, there exists no pure projection operator |^/')(^/'|, hence no positive operator 
V, which can be subtracted from PVgoj i leave them EWs again. So the EWs Wooi i are 
optimal. Similar argument proves the optimality of EWs WqI^^^^. 
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As for EWs Wi°^i_^i^, the state {ip) (if exises) ought to be of the form 

1^) = a++„k; +)\z; +)\z] -) + a^^Jz; +)\z; -)\z] +) 

[A. 2b) 

-)\z; +)\z; +) + a___\z; -)\z; -)\z; -). 

The same argument as above shows the impossibihty of existing such {tp). Therefore, the EWs 
j^j^ are also optimaL 

4.2 Optimality of EWs with conical FR 

The optimahty of EWs ^^^W^22i± has been proved in |36j, so we talk about the optimality of 
EWs ^^Wi22i±- Let us first find pure product states that the expectation value of "^'^'^^^122 i± 
over them vanishes. For this purpose, we consider a pure product state as follows 

W) = (g) (cos(|)|z;+) +exp(^¥,Jsin(|)|^;-)') (4.27) 
j=i V / 

and attempt to choose parameters 9. and if. such that 7''"[(^^'^Vr{2°2i-|-)|'^)('^|] = 0. By direct 
calculation, this trace is 



(4.28) 





= 1 ± COS 9-^ cos 


6*2 cos 6*3 H 


- sin 0^ sin 9^ sin 9^ 




X [cos COS cos(v32 + (■ 


-1)^+V3) 


+ sin ip sin (f-^ 


sin(</?2 + (- 


-1)V3)]- 


It is easy to see that the following four choices of parameters 


9. and o?. 


lead to zero 


the trace of ^^^W^i*2°2,o± 














^2 = ^3 = f ' 


n Stt 
^1 - — ' 
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^2 = ^3-- 
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~ 4' 
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^1 = ^3 = f 5 


n Stt 

"^2 ~ 2 ' 


Vi = ^, 


^2 = ^3-- 


TT 

~ 4' 
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^2 = ^3 = f ' 
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^1 - —■ 


, Vi = 


^2=^3-- 


TT 

~ 4 ' 




^1 = ^3 = h 


n Stt 

— 2 ' 


, Vi = -V", 


V2=Vz-- 


TT 

~ 4 ■ 
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For VI/^i22,o+; state (if exists) must be of the form f l4.25p and be orthogonal to the 



above four states, i.e., 



^+1' 


^) = iTlt- 


-«+++ 


+ 




^) = iTlt- 


~'^+++ 


— ia.^ 


{^,+ 




[-«++- 


1^ — ia 


{^4+ 




-«++- 


^ + ia^ 



. + 


-J] = 


= 0, 


. — ct_ 


-J] = 


= 0, 


+ a__ 


J] = 


0, 


— a 


J] = 


0. 


a. 


= a 


-+- 



The above system of four equations has trivial solution a^_^^ = = a_^_ = a _^ = 

provided that ip ^ ±|,±^. This proves the optimality of ^^^W^22o+ but ±|,±^ 

values of ip- Similarly, the optimality of '^^^W^22o- is proved for the same values of tp. 



5 Detection of p by EWs 

In this section, we consider the problem of detection of p by introduced EWs. 
5.1 Detection of EWs with polygonal FR 

First we begin with 16 EWs ^Wf" ^ ^ of (ES}. For these EWs we have 
TrCw[l,^,^p) = l + (-l)'-r333 + (-l)-'^r,,, + (-l)'3r,,, + (-l)"4r,,, + (-l)^^+^3+^4+V^^^. (5.29) 



It is seen that p is detectable by ^'^[J.^i^i^ if the parameters of p satisfy the following conditions 

b + c+- + -<±Ar^cosif., a + d + j + -<±Ar^cosif^, j = 1,2,3,4. (5.30) 
ad be 

For the 16 EWs ^W^^I^i^i^ of we have 

TrCw[l,^,^p) = l + (-l)Hr333 + (-l)^^r,,, + (-l)^3r^^^ + (_l)^4+V^^^ + (_l)«.+^3+^4^^^^. (5.3I) 



The detection condition imposes the following constraints on the parameters 

11 11 

b + c+ - + - < ±4r^smif., a + d + - + - < ±Ar^ simf., j = l,2,3,4. (5.32) 
ad J J be J J 
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5.2 Detection of EWs with conical FR 



The detection conditions obtained from 48 non- linear EWs F];-i.^^{p) of flXTSD together 
with 48 non-hnear EWs F^°i.^^{p) of (13141) are 

{a + l + h + \f<Aw , {a+\ + c+\f<Aw 

[b.66) 

{h+\ + d+\f <Aw , {c+\ + d + \f <Aw 
{a + \ + d + \f <4:W , {h+l + c+\f<Aw 
where w = u^^u^,u.^;v^,v^^v^ and 

= {r^ cosip^ ± cosf/Jg)^ + (r^ cosyjj =F ^4, cosip^)'^, 
= (r^ cos (f-^ ± Tg cos v^g)^ + cos -F cos (^4)^, 
Mg = (r^ cos ir^ cos (rg cos ¥?g ^r, cosy? J2, 

(5.34) 

f J = (rj sin ± Tg sin y9g)^ + (r^ sin yj^ =F sin (p^)^, 
^2 = (^1 sin V^i ± sin + (r^ sin =F sin V54)^, 
tig = (r^ sin ± sin yjj)^ + (rg sin p^ =]= sin p^y. 

5.3 Detection of EWs with cylindrical FR 

The detection conditions obtained from 36 non- linear EWs F^Ji.- ■ (p) of (13X7]) together 
with 36 non-hnear EWs F'^^i.-^-^{p) of flXTSj) are 



where 



< 16r^cos>^, < 16r^sinV,, i = 1,2,3; j = 1,2,3,4, (5.35) 



^ = (a + 5 + c + rf)(i + l + i + i), 

, = (a + 6 + i + i)(c + rf+i + i), (5.36) 



h = ia + c+l + ^)ib + d+^+' 



C I 



Unfortunately, as the following argument shows, the conditions (15.351) are not hold for p. We 
can write 

,a b. ,a c. ,a d. ,h c. ,b d. ,c d. 
4 + {- + -) + {- + -) + {- + -) + {- + -) + {- + -) + {- + - . 
a c a da c do d c 
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The two terms of each parenthesis are inverse of each other, so the value of each parenthesis 
is greater than or equal to 2 and hence 2;^ > 16, while in accord to fl5.35p < 16. Similar 
arguments show that z^jZ,^^ > 16, but in accord to (15.351) they are smaller than 16. 

5.4 Detection of EWs with spherical FR 

Finally, the detection conditions obtained from 18 non- linear EWs F^'-i.^{p) of ([321]) to- 
gether with 36 non- linear EWs ^'^"''^^5^^) of fl3:22|) are 

Zi<^u^ , z^<Av^, z,j = 1,2,3, 

where z^, and are defined as in (15.361) and (I5.34p . 



6 Comparison with other works 

If we put a = h = c = d=l, r^=r^ = l and ip^ = ip^ = 0, the detection conditions (15.331) give 

4 < 4 + (rg cosy^g — cosyj^)^. 

Hence, this case is detected by our EWs unless cosyjg = cosyj^. Further inspection shows 

that if in addition ip.^ = ip^ = 0, vr, then p is separable. So for the choice of parameters as 

a = h = c = d = 1, = = 1, ip^ = Lp^ = Q and p^ = p^ = 0, vr, the p is separable if and only 

if = r^; in agreement with Ref. [9]. 

For the case a = 1, < 6, c, ^ < 1, = 1, = and = r.^ = = 0, we have 

2(6 + c+i-3) 



2 + h + c + d+\ + \ + \- 



This trace attains its minimum value -0.3371 at 6 = c = ^ = 0.3798 and hence improves the 
result -0.1069 at 6 = c = i = 0.3460 of Ref. fT]. 



For the case a = 1, < ^, c? < 1, = 1, 1^9^ = and = = = 0, we have 

111 



2(i + H^-3) 
2 + 6 + c + rf+i + i + i- 



This trace attains its minimum value -0.3371 aX \ = - = d = 0.3460. 

c 
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7 2 2 (g) d Chessboard Density Matrices 

We generalize previous chessboard density matrices to 2 (8) 2 (8) d case and see that the PPT 
conditions are vahd. EWs forms remain the same with a few changes in notation. These 
methods can be apphed even for higher dimensions and for multi-qubits although the number 
of EWs and classification of them increases. Using some new algebraic notation for 2 (8) 2 (g) d 
case we can write 

1 d 
j=0 k=l 

+ z!f\OjP){l'ja\ + z'S\lja){0j(3\ + ^r|0j7)(lj7| 

ja ja jf 

+ dlj7)(0j7l + -^\lja){lja\ + (3) {Ij (3\ 

+ -Jt|1J7)(1J7| ) (7.37) 
here j = if j = 1 and vice versa and 

(3 = 0,...,d-l , <a < f3 <d-l , < < d - 1 

z_ = r_ expf t(p_ ) ,z_ = r_ expf —t(p_ ) 

For given a , /3 if < 1 for every j, ^, v then these type density matrices have positive partial 
transposes with respect to all subsystems, i.e., they are PPT states. All of previous witnesses 
classes including polygonal, conical, cylindrical and spherical become witnesses for this density 
matrices if we replace 

to { d X d identity matrix ) 
a to V2\\ 
a to ^A" 
a to E -E„„ 

on third partite of each terms of all of previous witnesses ( see appendix C ), as we do in 
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following subsections, where 

ay^ (3 = 0,...,d-l , 0<a<(3<d-l. 

7.1 Polygonal EWs 

With the notations as above, for polygonal case we have 32 (^^^^^) EWs. In analogy with 



([31]) the 16(^^) EWs are 

1 2 d 4. P (7.38) 

where E {0,1}^. The remaining 16(^^^^^) polygonal EWs can obtain by 

applying the phase-shift gate locally on the first qubit. The result is 

'W^:t3^4 = + i-^y^cr^^AE.. - EJ + V2{-iy^a^aX, 

7.2 Conical EWs 

We can expand the relevant density matrices in terms of Pauli and SU{N) operators for 
convenience ( see appendix C). In the following relations r.^^ are coefficients of relevant operator 
appearing in density matrices expansions, i.e. r.^^ is the coefficient of v^a.cr^.A^^, r..^ is the 
coefficient of \/2a.a.X . and r. .„ is the coefficient of a.a.(E^^ — E„„). 
The 96(^^^) conical EWs ( in analogy with fl3J3D and (EUD ) are 

(7.40) 

1 ± r,,,, - ^JX. + i-iyr.X + + Xy^X'- 

where k'j'l' is one of the triples 333, 330, 303, 033, and kjl is one of the triples 122, 212, 221. 



1 ±r 



(7.41) 
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here kjl is one of the triples 211, 121 and 112. Cyhndrical and spherical EWs for 2 ® 2 c? 
chessboard density matrices can construct with this procedure which are in full analogy with 
equations fl3l7l) , fl3A8D , (l3:2T|) , (13:221) . As result the number of EWs are 236(^^) 

7.3 2 2 3 Chessboard Density Matrices : An Example 

Now let us study density matrix for c/ = 3, Q; = 0,/9 = 2,7 = lin some details. In this case 
we can expand this density matrix in terms of Pauli and Gell-Mann operators , . . . , Ag ( 
see appendix D ), and all of previous witnesses including polygonal, conical, cylindrical and 
spherical are valid if we replace 
Jj to /3 ( 3 X 3 identity matrix ) 

to y2A+ = A, 

to V2X;^ = A, 
a, to ^00-^22 = |(A3 + V3A3) 

on the third partite of each terms of all of previous witnesses. For example, using above 
prescription, polygonal witness in (13. Sp can be written as 

+ (-l)^3a a A, + i-lY^a a A, + (-l)*2+»3+»4+V a A, 

'V / X y 5 \ / y X 5 \ / y y 4: 

By similar substitution, all of 236 EWs can be constructed. The detection ratio ( the ratio of 
entangled density matrices detected by all our EWs to all randomly selected density matrices 
), is listed in table 2. 



8 Numerical analysis of entanglement property of p 

In this section we deal with some numerical analysis regarding detection ability of introduced 
EWs for 2® 2® 2 and 2® 2® 3 chessboard density matrices. Numerical calculation is done on 
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random set of relevant PPT chessboard density matrices. Those density matrices detected by 

EWs are counted and then the ratio is calculated. The percent of the volume of phase space 
that can be detected by introduced EWs is as listed in the table 1. 



EWs 


percent of detection 


EWs 


percent of detection 


Polygonal 


28.3 


Not polygonal but conical 


0.44 


Conical 


18.3 


Not polygonal but spherical 


0.0275 


Spherical 


0.047 


Polygonal and spherical 


0.0176 


All EWs 


28.62 


Conical and spherical 


0.031 



Table 1: The percent of detection for introduced EWs. "Not polygonal but conical" means the 
percent of the three-qubit PPT density matrices p that the polygonal EWs can not detect but conical 
ones can detect. 



EWs 


percent of detection 


AU 236 EW's 


R ± a = 85.45 ±3.336 



Table 2: The percent of detection for introduced 2 (g) 2 (g) 3 EW's. R indicates mean ratios and a is 
standard deviation 

9 Conclusion 

In this paper, we have considered a class of three-partite PPT Chessboard density matrices 
and via an exact convex optimization method, have constructed various linear and non-linear 
EWs detecting them. The operators participating in constructing the EWs have been chosen 
such that the geometrical shape of the feasible region have been obtained exactly. The EWs 
have been classified according to the geometrical shape of relevant feasible regions. When 
feasible region was not a polygon, non-linear EWs were obtained. The optimahty of EWs 
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with polygonal and conical feasible region have been shown. The introduced EWs were all 
non-decomposable, since they were able to detect PPT entangled states. Event hough, we 
have mainly discussed these methods for 2 (g) 2 (8) 2 and 2 (g) 2 (8) o? chessboard density matrices, 
but they are general and one can apply them for di <^ d2 <^ ds via some minor changes in 
notation and calculations. It was shown that the detection ability of introduced EWs is often 
comparable with one of EWs introduced elsewhere. In some cases, the detection abihty of EWs 
introduced here is better. Finally the prescription of this work is applicable for multi-partite 
PPT Chessboard density matrices which is under investigation. 
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Appendix A 

Proving the inequalities: 

In the foUowing proofs, we use the abbreviations 

Tr(af^ \a){a\) = 

Tr(af) \/3){(3\) = b^ (A-i) 
rr(af |7)(7l) = c,. 

Since + + = 1 and also the similar relations hold for 6.'s and C-'s, so the points a, b, c 
lie on a unit sphere and we can parameterize their coordinates by using spherical coordinates 
6 and as follows 

= sin6'j cosyj^, = sin 6'^ sin*/?^, = cos^^^ 
6^ = sin cos , = sin sin (/^^ , b^ = cos 
= sin 6*3 cos , C2 = sin 9^ sin , 03 = cos 9^ . 

The proof of SSTTh : 

To prove this equality, we note that 

Pi = a^b^c^ = cos 9-^ cos 9^ cos 6*3 

P2 = «i(&iCi ± 62C2) = sin6'j sin 6^2 sin 6*3 cosip-^ '^os{(p.^ =F ^2) 
P3 = a^ip^c^ -F &2'^i) = sin 6'^ sin 6^2 sin 6^3 siny^^ sin((y93 =f ^92) 



whence 



p2 p2 

^ + — ^ — = sin^ sin^ 6^2 sin^ 9^ 



cos^ sin^ 
Taking derivative with respect to ip^ we obtain 

P2 , P3 



Above two equations yield 



cos^ V?^ sin^ yj^ 



. 2 ±P3(P2±P3 

sm 



sin 9^ sin ^^2 sin 9^ 
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2 P2{P2±P3 
COS if^ 



sin 9-^ sin 9^ sin 

Noting that sin^ ip-^ + cos^ = 1, we get 



A + -P3 = ± sin sin 9.^ sin ^3 , P2 — -P3 = =t sin ^1 sin sin 9^ 



Eliminating between P2 ± -P3 and Pi leads to 



sin^ sin^ 6^3 cos^ 9^ cos^ 6^3 



Taking derivative with respect to 9^ and by similar argument as above, we get 

Pi P2±P3 Pi P2±P3 

' ^ ±1 , — ^— ^ = ±1. 



COS ^3 sin ^3 cos6'3 sin 6^3 

Finally, taking derivative with respect to 9^ and using the identity sin^ 9^ + cos^ ^3 = 1 gives 

Pl^ + (P2±P3)i = l 

But, as the Fig. 1 shows, this is a concave curve. Since the mixed separable states are convex 
combinations of pure product states, the boundaries of FR are the planes of (13.71) . 
The proof of (IXTOb : 

The proofs are similar, so we give the proof for the case Qf" = O333 and kjl = 122. We note 
that 

Pi = aj)^c^ = cos^^ COS02 cos ^3, 
P2 = ai(&iCj ± &2C2) = sin^^j sin^^2 sin 6*3 costy?^ cos^Lp^ =1= ^93), 
P3 = a^{bj^c^ ± 62C1) = sin6'j sin 6*3 sin 6*3 sin ip-^ sin((^2 ^fa)- 
By eliminating 9-^ and (p-^ , we get 

, 1 f PI , PI 



cos^ 6^2 cos^ 9^ sin 9^ sin 9^ \cos2(v92 =F V's) sin (932 
Now we put V^2 = V^3 = f = ^ ^nd ^93 = | to obtain 

p2 p2 _|_ p2 

-'1 , 2 "T 3 1 



cos^ 6*2 cos^ 6*3 sin^ 9^ sin^ 6'3 
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Derivation with respect to 6^ leads to 



Pi ^ (Pi + P? 



cos^ cos sin 6^ sin 9^ 



Above two equations yield 



sm \ cos sm 6^ I cos \ cos 9^ sm fe^g 

From sin^ 9^ + cos^ = 1, we have 

Pi ^{Pi + PDi^^^ Pi + 



cos ^3 sin ^3 cos ^3 sin ^3 

Finally, taking derivative with respect to 9^ and using the identity sin^ 9^ + cos^ ^3 = 1 gives 

But, as the Fig. 1 shows, this is a concave curve in terms of variables Pi and {P2 + P^)^ ■ 
Since the mixed separable states are convex combinations of pure product states, the relations 
between these two variables are given by the lines 

Pi + (P| + P|)5 = ±1, Pi - (P| + P^)-2 = ±1 

So the relations between Pi, P2, and P3 are as in (13.101) . 

If we take Qi° = O330, the proof of (I3.10p proceeds as follows. We note that in this case 

Pi = a^h.^ = cos cos 

By eliminating and ip-^ , we get 

P? , 1 f Pi , Pi 



+ 



008^9 2 sin sin ^^3 \cos'^{(p2 T ^3) sin (^92 ±7^3) 
Now we put 6^3 = I and '/'a = V^a = f V^2 = t '/'a = f obtain 

P? ^Pi + Pi_, 



cos^ 9^ sin^ 9^ 
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Taking derivative with respect to 6^ leads to 

cos^ sin^ 

Above two equations yield 

sin^ = ±(P| + P|)i(Pi ± (P| + P|)^), cos^ 9, = Pi(Pi ± (P| + P|)^) 
Finally, the (13.101) follows from the identity sin^ 9^ + cos^ 9^ = I- 
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Appendix B 

The coefficients of Pauli operators appearing in p: 



^11 


2 

= n< 


cos </?j + cos (/?2 + COS (/?3 + cos (^^ j 


^300 


1 

n 


[a 


1 U 1 


^ 1 ^ 
c + a — 


1 1 

a b 


1 1 
c d 


'"112 


_ 2 
n 


(r^ sin - sin ip^ + sin ip^ - sin (/?^) 


^030 


1 
n 


[a 


+ — 


c — a — 


1 1 I 
a 6 


1 I 1 
c d 


T 

' 121 


2 

n 


(rj sin (^^ + sin (p^ — sin ip^ — r\ sin 994) 


^003 


n 


a 


-b + 


c — d — 


1 1 

- + - - 

a ^ b 


- + - 


' 211 


n V 


— sin — sin ip^ — sin ip^ — sin (/j^) 


r — 

' 330 


n 


[a 


+ b- 


c-d + 


b 


1 _ 1 
c d 


r = 
' 122 


= -(- 


-r^ cos + cos 'p.^ + r3 cos (p.^ — r\ cos (/:'4) 


r = 

303 


n 


a 


-b + 


c-d + 


a b ^ 


1 _ 1 

c d 


' 212 


= H 
n ^ 


Tj COS (^j — COS + COS ip^ — COS f^^) 


' 033 


n 


[a 


-b- 


c + d + 


1 _ 1 _ 

a b 


d 


' 221 


= 

n ' 


COS </7j + COS (/?2 — Tg COS ip^ — COS (/J^) 


r — 
' 333 


n 


[a 


-b- 


c + d- 


- + - + 


1 _ 1 
c d 


r 

222 


_ 2 
n 


{r\ sin (/)^ — sin ^p^ — sin (^3 + sin 994) 
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Appendix C 



Every d-dimensional square matrix could be written in terms of square matrices E.^, which 
show the value 1 at the position and zeros elsewhere. Now one can define Hermitian 
traceless basis for d-dimensional matrices as follows ( see |3i 
The off-diagonal basis are 



A 



a/3 



jV2 



and diagonal basis are 



0/ 



V3 



d{d-l) 



( 



Vo 0/ 

In order to generalize the witnesses, we must write E^^ in terms of { d x d identity matrix 
) and A„'s. Some calculation shows that 



E-- = E- . . - . , -\- \ I — 7 r A. — 



2(z + 1) 



A. , , 0<i< d-2 

2{i + l) '-^ ' - - 



recursion relation ) and 



E = -I 



d-l , 
~2d' 



Proving the inequalities for 2 2 ® 
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The proof is almost the as explained in appendix A. We use the abbreviations 

rr(<7f) \a){a\)^a, 
Tr{af^ \(3){(3\)^h, 
Tr{V2Xl |0(el) = q (A-ii) 

where 

10 



2 + + r2 , 

d-l \ r^_^e 



We have 



12 3 '12 3 



and 



/• 2 _^ 2\2 
2 _|_ 2 _|_ 2 g ' 13' 

•^1 + "^2 + "^3 " TIS^T 772"^ ~ ^ 
if we set g = 1 without loss of generality, then the points a,b,c lie on a unit sphere and we can 
parameterize their coordinates by using spherical coordinates 9 and (f as follows 

— sin 9^ cos </7^ , — ^1 V'n ^3 — ^^^^ ^1 
6^ = sin ^2 cos </?2 , &2 = sin ^2 sin (p^, b^ — cos ^2 
= sin 6*3 cos ip^, c^= sin 6*3 sin 993 , C3 = cos 6*3 . 
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Appendix D 

The Gell-Mann Matrices 



The analog of the Pauh matrices for SU(3) are GeU-Mann matrices defined as: 
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Figure Captions 

Figure-1: The boundaries of feasible region for pure product states (dotted curve) and 
mixed separable states (line) for EWs of relation (13.81) . 
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